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ABSTRACT 

An  analysis  of  steady  state  operation  of  a  vertically 
oriented,  variable  conductance  heat  pipe  is  presented.   The 
effects  of  binary  mass  diffusion,  axial  pipe  wall  conduction, 
and  gravitational  effects,  caused  by  the  difference  in 
molecular  weight  of  the  non-condensible  gas  and  the  working 
fluid  are  incorporated.   Analytical  expressions  for  the  con- 
servation of  mass,  momentum,  and  energy  are  combined  along 
with  equations  of  state  to  describe  steady  state  operation 
of  the  heat  pipe.   These  expressions  are  combined  to  form 
a  system  of  three  differential  equations  with  three  unknowns; 
working  fluid  vapor  velocity,  mass  fraction,  and  temperature. 
These  equations  are  nondimensionalized  and  an  iteration  scheme 
for  numerical  solution  of  the  equations  is  presented. 
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I.   INTRODUCTION 

A.   BACKGROUND 

A  heat  pipe  is  a  closed,  evacuated  chamber  that  has  the 
inside  walls  lined  with  a  wick  structure,  saturated  with  a 
volatile  working  fluid.   The  steady  state  operation  involves 
vapor  heat  transfer  and  capillary  action.   A  comprehensive 
coverage  of  the  theory  and  operation  of  heat  pipes  can  be 
found  in  Dunn  and  Reay  [Ref.  1].   Briefly,  the  operation  of 
a  vertically  oriented  heat  pipe  will  be  discussed.   The  heat 
pipe  has  three  sections:   an  evaporator  section,  an  adiabatic 
section,  and  a  condenser  section.   The  heat  pipe  with  the 
three  sections  marked  is  shown  in  Figure  1.1.   The  lower  evapora- 
tor end  of  the  heat  pipe  is  heated,  causing  the  working  liquid 
in  the  wick  to  vaporize.   This  phase  change  causes  the  working 
fluid  to  absorb  latent  heat  and  the  vapor  pressure  increases 
in  the  evaporator  end.   The  working  fluid  is  drawn  away  from 
the  wick,  causing  the  liquid  pressure  to  decrease.   The  vapor 
moves  to  the  cooler  condenser  where  it  is  condensed  and  the 
working  fluid  gives  up  its  latent  heat.   The  condensate  is 
deposited  on  the  wick  which  causes  the  liquid  level  to  in- 
crease thus  increasing  the  liquid  pressure.   Capillary  forces 
take  the  higher  pressure  liquid  from  the  condenser  to  the  lower 
pressure  liquid  end  of  the  evaporator.   The  capillary  forces 
are  assisted  by  gravitational  forces  when  the  heat  pipe  is 
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Condenser 


'    Adiabatic  Section 


Figure  1.1.   Heat  Pipe 
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vertically  oriented  as  shown  in  Figure  1.1.   It  is  important  to 
note  that  for  heat  pipes  the  evaporator  position  is  not  re- 
stricted, but  since  the  experimental  data  and  numerical  model 
use  a  vertically  oriented  heat  pipe,  this  orientation  will  be 
used  throughout  this  thesis. 

The  great  advantage  of  heat  pipes  is  that  the  vaporization 
and  condensation  process  allow  for  high  heat  transfer  rates 
with  small  temperature  differences.   Two  other  advantages  are 
that  the  basic  heat  pipe  has  no  moving  parts  and  can  be  used 
in  many  geometrical  configurations.   The  heat  pipe  has  limi- 
tations and  again,  a  comprehensive  presentation  of  these  limi- 
tations can  be  found  in  Dunn  and  Reay  [Ref.  1].   Briefly,  the 
five  heat  pipe  limitations  are  as  follows: 

(a)  The  viscous  limit  is  basically  the  low  temperature  limit, 
This  would  be  the  minimum  temperature  for  the  heat 

pipe  to  begin  operation.   Generally,  the  viscous  limit 
is  only  important  during  start  up  of  some  heat  pipes, 
i.e.,  liquid  metal  pipes. 

(b)  The  sonic  limit  describes  the  limit  of  the  vapor  flow 
velocity  at  the  heat  pipe  evaporator  exit.   The  re- 
sulting choked  flow  condition  limits  the  axial  heat 
flux  capacity  of  the  heat  pipe. 

(c)  The  entrainment  limitation  is  due  to  the  removing  of 
the  liquid  before  it  can  return  to  the  evaporator  by 
vapor  flowing  at  high  velocity. 

(d)  The  wicking  or  capillary  limit  reflects  the  fact  that 
each  combination  of  wick  and  fluid  has  a  maximum 
capillary  head.   As  stated  earlier,  the  capillary  head 
must  increase  with  the  liquid  and  vapor  pressure  drops 
and  can  do  so  only  up  to  the  maximum  capillary  head. 

(e)  The  last  limit  is  the  boiling  limit.   If  temperatures 
in  the  evaporator  get  too  high  the  working  fluid  will 
begin  to  boil  within  the  wick.   This  will  disrupt  liquid 
flow  to  the  evaporator,  causing  the  heat  pipe  to  stop 
working . 
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These  five  heat  pipe  limitations  vary  in  importance,  depending 
on  the  temperature  and  corresponding  axial  heat  flux  operating 
point.   Since  at  different  temperatures  one  can  encounter 
different  heat  pipe  controlling  limitations,  it  is  important 
to  have  a  basic  understanding  of  all  five  limitations. 

B.   VARIABLE  CONDUCTANCE  HEAT  PIPE 
1 .   Basic  Operation 

Variable  conductance  heat  pipes  are  described  in  detail 
in  Marcus  [Ref.  2].   As  described  in  the  background  section,  a 
change  in  the  heat  load  will  automatically  change  the  operating 
temperature  and  pressure  of  the  heat  pipe,  given  a  constant  sink 
temperature.   This  is  an  undesirable  quality  for  many  situations 
where  it  is  desired  to  have  a  stable  source  temperature  under 
conditions  of  varying  heat  load  or  sink  temperatures.   An 
example  of  this  case  would  be  the  cooling  of  sensitive  elec- 
tronic equipment.   Basically,  a  noncondensible  gas  is  introduced 
into  the  heat  pipe  along  with  the  working  fluid,  making  the 
heat  pipe  become  a  passively  controlled  device.   The  nonconden- 
sible gas  occupies  a  certain  portion  of  the  condenser  section, 
depending  on  operating  conditions.   Increasing  operating  tempera- 
ture in  turn  increases  the  vapor  pressure  of  the  working  fluid 
which  compresses  the  noncondensible  gas  into  a  smaller  volume. 
This  provides  a  greater  active  condenser  area,  thus  allowing 
a  m.ore  stable  source  temperature  than  a  conventional  heat  pipe. 
When  the  operating  temperature  decreases,  working  fluid  vapor 
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pressure  falls  and  the  noncondensible  gas  expands,  blocking  a 
larger  portion  of  the  condenser.   The  result  is  a  passively 
controlled  variable  condenser  area,  which  in  turn  makes  a 
more  stable  source  temperature. 

2 .   Variable  Conductance  Heat  Pipe  Theory 

Marcus  [Ref.  2]  describes  a  flat  front  theory  and  a 
'diffuse  front  theory'  for  variable  conductance  heat  pipes. 
Flat  front  theory  assumes  the  noncondensible  gas  is  in  the 
form  of  a  plug  with  sharply  defined  vapor-gas  interface  which 
is  perpendicular  to  the  axis  of  the  pipe.   Since  this  model 
neglected  the  axial  conduction  in  the  heat  pipe  wall  and  diffu- 
sion of  vapor  into  the  stagnant  gas  plug,  experimental  results 
[Ref.  2]  have  shown  this  model  to  be  inaccurate.   The  early 
model  of  flat  front  theory  provided  initial  insight  into  the 
operating  characteristics  of  gas  loaded  heat  pipes.   With  the 
incorporation  of  binary  diffusion  between  the  working  fluid 
vapor  and  the  noncondensible  gas  and  axial  conduction  in  the 
pipe  wall,  Marcus'  [Ref.  2]  'diffuse  front  theory'  establishes 
a  model  of  gas  loaded  heat  pipes  which  provides  better  simi- 
larity between  computer  program  and  laboratory  results. 

Gravitational  effects  on  the  vapor-gas  interface  are  not 
included  in  flat  front  theory  or  the  'diffuse  front  theory.' 
Gravitational  effects  can  greatly  distort  the  axial  temperature 
profile  when  the  working  fluid  and  noncondensible  gas  have 
significantly  different  molecular  weights.   This  fact  is  demon- 
strated by  experimental  results  reported  by  Kelleher  [Ref.  3]. 
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Batts  [Ref.  4]  conducted  variable  conductance  heat  pipe  experi- 
ments using  working  f luid/noncondensible  gas  combinations  of 
methanol/helium  and  raethanol/krypton .   Kane  [Ref.  5]  developed 
an  analytical  model  of  variable  conductance  heat  pipes  which 
includes  the  effects  of  binary  mass  diffusion,  axial  pipe-wall 
conduction,  and  gravitational  effects  which  exist  at  the  vapor- 
gas  interface  region  as  shown  by  Batts  [Ref.  4]  experiments. 
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II.   FORMULATION  OF  AN  ANALYTICAL  MODEL 

A.  PROBLEM  STATEMENT 

The  purpose  of  this  chapter  is  to  formulate  the  equations 
describing  the  steady  state  operation  of  a  vertically  oriented, 
gas  loaded,  variable  conductance  heat  pipe.   As  in  Figure  1.1, 
the  condenser  is  above  the  evaporator  for  this  orientation. 
This  analysis  includes  the  effects  of  binary  mass  diffusion, 
axial  pipe  wall  conduction,  and  gravitational  effects  due  to 
the  difference  in  m.olecular  weight  between  the  noncondensible 
gas  and  the  working  fluid. 

B.  ASSUMPTIONS 

The  following  assumptions  are  used  to  develop  the  system 
of  equations: 

1.  Heat,  mass  and  momentum  transport  processes  are  treated 
as  one-dimensional  in  the  axial  direction  within  the 
vapor  space  of  a  vertically  oriented  heat  pipe. 
Kelleher  [Ref.  3]  shows  that  in  any  other  orientation, 
gravitational  effects  cause  highly  three-dimensional 
spatial  variations  when  the  molecular  weights  of  the 
working  fluid  and  noncondensible  gas  differ  substantially 

2.  The  noncondensible  gas  acts  as  a  perfect  gas.   The 
mixture  of  noncondensible  gas  and  working  fluid  are 
treated  as  thermodynamically  independent  substances. 

3.  Because  of  direct  contact  between  the  vapor  and  the 
liquid  saturated  wick,  the  working  fluid  vapor  is  at 
saturation  conditions  at  all  times. 

4.  The  thermal  resistance  of  the  wick-wall  structure,  in 
the  radial  direction,  is  negligible  in  comparison  with 
the  thermal  resistance  due  to  convection  outside  the 
pipe. 
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C.   ANALYSIS 

Kane  [Ref.  5]  presents  a  proposed  solution  for  this  prob- 
lem.  The  following  is  a  significant  change  to  his  proposed 
solution  with  modifications  being  noted  as  they  occur  in 
this  section.   Variables  with  the  subscript  'a'  apply  to  the 
working  fluid  and  variables  with  the  subscript  'b'  apply  to 
the  noncondensible  gas. 

Utilizing  the  assumption  that  the  noncondensible  gas  is 
a  perfect  gas,  the  equation  for  noncondensible  gas  pressure 
is  : 


Pb   =   pb  i^T  <"•!> 


The  principles  of  mixtures  of  thermodynamically  independent 
substances,  which  are  found  in  Reynolds  and  Perkins  [Ref.  6], 
are  used  to  define  the  total  density  and  the  total  pressure. 


Pa  -^  Pb  (^^-2) 


P   =   P^  +  P,  (II. 3) 

a    b 


The  definition  of  mass  fraction  is  used  to  form  the  equation 
for  the  working  fluid  vapor  mass  fraction. 


Pa 
CO    =   —  (II. 4) 

a      p 
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The  solution,  utilizing  the  conservation  of  mass,  is 
based  on  a  differential  element  representing  a  section  of 
the  heat  pipe  condenser,  Figure  2.1.   The  conservation  of 
mass  of  species  'a'  for  the  differential  element  yields: 


pVA    =   (pV   +:r— (pV)dx)A^+mPdx. 
^a  a  c      ^a  a   dx   a  a     c     a  v 


This  equation  sets  the  mass  of  species  'a'  into  the  control 
volume  equal  to  the  mass  of  species  'a'  out,  including  a 
term  to  account  for  the  removal  of  species  'a'  by  condensation, 
The  process  of  condensation  is  accounted  for  by  multiplying 
the  condensation  mass  flux  times  the  area,  defined  by  the 
vapor  space  perimeter  times  the  length,  dx ,  of  the  differen- 
tial element.   After  basic  algebra,  the  equation  for  conser- 
vation of  mass  of  species  'a'  becomes: 


0   =   mP   +A   ^vpV' 
a  V    c  dx   a  a 


By  using  Fick's  law  of  binary  diffusion  [Ref.7],  and  solving 

for  p  V  ,  the  resulting  formula  is: 
a  a 

doj 

p  V    =   00  pV  -  pD  —5 —  . 
•^a  a      a^     ^   dx 


The  final  form  for  the  conservation  of  mass  for  species  'a' 
becomes : 


P     ,  ,      do) 

m   ^  +  §-(lo  pV)  -  §-(pD  -v-^)   =   0  .  (II. 5) 

a  A     dx   a^      dx  '^   dx 
c 
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since  species  'b'  is  a  noncondensible  gas,  the  condensation 
mass  flux  of  species  'b'  equals  zero.   The  resulting  equation 
for  the  conservation  of  mass  for  species  'b'  is: 


,  -,      doo, 

^(co,  pV)  -  §-(pD  -T^)   =   0 
dx   b^      dx  '^   dx 


The  mass  continuity  of  the  mixture  is  derived  by  adding 
the  conservation  of  mass  of  species  'a'  with  the  conservation 
of  mass  of  species  'b'. 


P     -,  ,      d(jo     dco, 

m   ^--  +  3— (pV(u3  +03,  ))  -  :g— (pD(— 5 —  +  —5 — )   =   0  . 
a  A     dx      a    b      dx^   dx     dx 
c 

Because  go  +co,  is  equal  to  one,  and  :5— (co  +  go,  )  =  35— (1)  =  0, 
a    b      ^  '      dx   a    b     dx        ' 

the  final  result  of  the  mass  continuity  of  the  mixture  is: 


c 


Figure  2.2  is  a  differential  element  of  the  vapor  space  in 
the  condenser  and  is  used  to  derive  the  conservation  of  linear 
momentum  equation  for  the  mixture  of  noncondensible  gas  and 
working  fluid  vapor.   The  conservation  of  linear  momentum  for 
the  mixture  of  noncondensible  gas  and  working  fluid  vapor  is 
derived  by  setting  the  summation  of  external  forces  acting  on 
the  control  volume  in  the  x  direction  equal  to  the  net  efflux 
of  linear  momentum. 


If    =  /      V(pV  dA) 
CS,x 
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The  summation  of  external  forces  acting  in  the  x  direction  are 


F         =       pgAdx-TPdx-dpA 
X  ^         c  V  ^      c 


These  terms  represent  the  force  due  to  gravity  acting  on  the 
control  volume  minus  the  forces  due  to  wall  shear  stress  and 
pressure  acting  on  the  control  volume.   The  wall  shear  stress 
can  be  calculated  by: 


X   =   >|v2 


where  A  is  a  dimensionless  friction  factor.   If  developed 

laminar  flow  is  assumed  the  Darcy  friction  factor,  A  =  :^— , 

Re 

could  be  used  with  the  Reynolds  number  defined  as 

p  VD 

V 

Re   = 

y 
The  equation  for  wall  shear  stress  then  becomes: 


V 


where  y  is  the  viscosity  of  methanol  vapor,  V  is  the  velocity 
of  the  vapor  and  D   is  the  diameter  of  the  vapor  space.   The 
terms  to  represent  the  net  efflux  of  linear  momentum  are: 


j      V(pvdA)   =   V  m  P  dx  +  (pv  +d(pv)  (v+dv))A^  -  pV  A 
CS,x  a  a  V  c       c 

=   V  m  P  dx  +  A  pvdv  +  A  vd(pv)  . 
a  a  V       c         c         *^ 
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By  setting  the  suirmation  of  external  forces  in  the  x  direction 
equal  to  the  net  efflux  of  linear  momentum  and  combining  like 
terms,  the  result  is: 


dp  d  -   2,     V  _  a     V  Syv  ,^^    _,, 

^  =  pg  -  ^(Pv  )  -  _v—  -  ^  ^  (II. 7) 

c    a    c   V 


The  above  equation  represents  conservation  of  linear  momentum 
for  the  mixture  of  noncondensible  gas  and  working  fluid  vapor. 
Kane  [Ref.  5]  assumes  zero  velocity  for  the  noncondensible  gas 
and  derives  an  equation  for  conservation  of  linear  momentum 
of  species  'b.'   The  species  'b'  conservation  of  linear  momen- 
tum equation  is  not  necessary  to  derive  the  model  in  this 
analysis.   Since  the  conservation  of  linear  momentum  for  the 
mixture  of  the  noncondensible  gas  and  working  fluid  vapor  is 
the  more  consistent  approach.  Equation  (II.  7)  is  used  through- 
out this  thesis. 

Figure  2.3  is  a  differential  element  of  the  condenser  wall 
and  represents  the  wick  and  pipe  wall  combination  as  one  ele- 
ment.  The  wick-wall  combination  is  treated  like  a  fin  and  is 
used  to  derive  the  equation  for  the  conservation  of  energy. 
Heat  transfer  is  represented  by  the  latent  heat  deposited  on 
the  'fin'  by  condensation  of  the  working  fluid  vapor  plus  the 
heat  transfer  from  Fourier's  law  of  conduction  entering  the 
'fin. • 


Q.    =   m  h^  P  dx  +  Q 
in       a   fg  V        c 
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m      h^     P    dx 
a      fg     V 


h  P     (T  -T    )dx 
s  °° 


Q     +dQ      = 
c  c 


,.^       dT    ^    d     ,    -.   ^     dT>  , 

•KA      3—  +   :5— (-K  A    :5— )dx 

p   dx        dx  pdx 


Figure  2.3.   Differential  Element  for  Energy  Balance 
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Heat  transfer  out  is  represented  by  Fourier's  law  of  conduction 
leaving  the  'fin'  plus  the  energy  which  is  convected  away  from 
the  fin. 


Q^.   =   Q   +dQ   +hP(T-T)dx 
out       c      c       s      °° 


The  equation  for  conservation  of  energy  is  derived  by  setting 
heat  transfer  in  equal  to  heat  transfer  out  of  the  'fin.' 


^2^  hP  h^  P   . 

^  -  ^(T  -T  )  +  42_v:m    =   0             (II. 8) 

-,  2  kA       oo  kA     a 

dx  p  p 


The  experimental  data  from  Batts  [Ref.  4]  is  used  to 
correlate  this  analysis.   The  working  fluid  is  methanol, 
CH-,OH.   The  working  fluid  vapor  in  the  heat  pipe  is  in  direct 
contact  with  the  working  fluid  liquid  in  the  wick  and  is 
therefore  considered  to  be  a  saturated  vapor  at  the  local 
temperature  of  the  wick-wall  combination.   Eubank  [Ref.  8] 
tabulates  data  for  the  saturation  properties  of  methanol. 
These  tables  are  used  to  generate  functions  of  local  temperature 
for  partial  pressure,  partial  density  and  heat  of  vaporization 
of  saturated  methanol  vapor. 

(a)   The  partial  pressure  of  saturated  methanol  vapor  is  a 
function,  Fa,  of  local  temperature. 


P    =   Fa(T)  (II. 9) 

a 
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(b)   The  partial  density  of  saturated  working  fluid  vapor 
is  a  function,  Sa,  of  local  temperature. 


p^   =   Sa(T)  (11.10) 

a 


(c)   The  working  fluid  heat  of  vaporization  is  a  function. 


H^  ,  of  local  temperature 


h^    =   H^  (T)  (11.11) 

f  g      f  g 


In  summary,  the  eleven  equations  which  are  stated  in  this 
chapter  form  the  foundation  for  the  analysis  of  a  vertically 
oriented,  variable  conductance  heat  pipe.   The  remaining 
chapters  derive  a  numerical  model  and  iteration  scheme  based 
on  the  eleven  equations  from  this  chapter. 
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III.   NUMERICAL  MODEL  FORMULATION 

A.   VARIABLE  ELIMINATION 

The  system  of  eleven  equations  can  ultimately  be  reduced 
to  a  system  of  three  equations  by  eliminating  variables. 
Equation  (II. 9)  is  eliminated  by  using  the  equation  of  state, 

Fa(T),  for  P  .   The  partial  density  of  the  saturated  working 

a 

fluid  vapor,  p   (Equation  11.10)  is  eliminated  by  using  Sa(T). 
H^  (T) ,  Equation  (11.11),  is  used  to  eliminate  h^    .   The  non- 
condensible  gas  pressure,  P,  ,  is  eliminated  by  using  the  perfect 
gas  equation  of  state  (Equation  II. 1).   The  treatment  of  working 
fluid  vapor  and  noncondensible  gas  as  thermodynamically  inde- 
pendent substances  gives  total  pressure  as  the  sum  of  partial 
pressures  and  total  density  as  the  sum  of  partial  densities. 
Total  pressure  is  eliminated  by  using  the  sum  of  ^partial  pres- 
sures. Equation  (II. 3).   The  noncondensible  gas  density  is 
eliminated  by  solving  Equation  (II. 2)  for  p,  and  replacing  p, 
with  p  - p  .   The  definition  for  mass  fraction  of  the  working 
fluid  vapor.  Equation  (II. 4),  is  used  to  replace  total  density. 

After  making  the  preceding  seven  substitutions,  four  equa- 
tions with  four  unknowns  are  left.   The  four  remaining  depen- 
dent variables  which  are  functions  of  axial  position  are: 

2 

(a)  m    ,    working    fluid   condensation    flux    [Ibm/ft      sec] 

a 

(b)  V,  velocity  of  the  mixture  [ft/sec] 

(c)  CO  ,  working  fluid  mass  fraction  [0  <  to   <  1] 

a  a 

(d)  T,  temperature  [°R]. 
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The  four  remaining  equations  are 


c  a 


"^'      "a  A^-lr'Sa  V)    -§3^(^^)      =      0  (III. 2) 

c  a 


a  b 


=       (S£  g)         d      Sa   ^2  ^^:^      ^    8y    ^  3 

00       ^  dx    CO  A         CO         A       D 

a  a  c         a        c      V 


«3)       f^-    ifr  h(T-TJ     +    ^^^       =       0  (III. 4) 

dx  p  p 


Furthermore,  Equation  (III.l)  can  be  eliminated  by  solving  this 

equation  for  working  fluid  condensation  flux,  m  ,  and  substi- 

a 

tuting  that  into  the  three  remaining  equations. 

The  resulting  three  equations  with  dependent  variables  of 
velocity,  working  fluid  mass  fraction,  and  temperature  are: 


(a)       ~(—  V)     +    5-(Sa    V)     -    g-(^D-^)       =       0  (III. 5) 

dx   CO         '         dx  dx   CO  dx 

a  a 


,,  X       dFa    ,    d     ,^     ,   1  ,  ,      R    m\ 

(b)  -5 +     3—     Sa 1)       rr—     T 

dx         dx           CO  M, 

a  b 

=       (Sa      )  d      Sa      2      ^^d      Sa^         V^^ 

CO      ^  dx    CO                    CO      dx    CO                 AD 

a  a                  a             a                c   v 
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ax       p  pa 


This  set  of  three  equations  and  three  unknowns  describe  the 
operation  of  a  vertically  oriented, gas  loaded, variable  conduc- 
tance heat  pipe. 

B.   BOUNDARY  CONDITIONS 

These  equations  describe  the  thermodynamics  of  the  heat 
pipe  condenser.   The  independent  variable,  the  distance  along 
the  heat  pipe,  is  measured  with  the  origin  at  the  condenser 
end  farthest  from  the  evaporator.   The  total  condenser  length 
is  defined  as  L.   Based  on  the  system  of  three  equations,  two 
boundary  conditions  are  required  for  temperature  and  mass  frac- 
tion and  one  boundary  condition  is  required  for  velocity. 

From  experimental  data,  the  temperature  at  the  far  end  of 
the  condenser  and  in  the  evaporator  are  known.   The  boundary 
conditions  for  the  dependent  variable  of  temperature  are: 


when  X  =  0,  T  =  T  ,  where  T   equals  the  temperature  at  the 

far  end  of  the  condenser. 

when  X  =  L,  T  =  T  ,  where  T   equals  the  temperature  at 

the  end  of  the  condenser  closest  to 
the  evaporator. 


The  highest  velocity  of  the  working  fluid  vapor  is  reached 
in  the  evaporator  exit.   This  velocity  can  be  calculated  from 
the  known  heat  input  to  the  heat  pipe.   The  vapor  velocity  in 
the  evaporator,  V  ,  is  equal  to  the  heat  transfer  rate  into 
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the  evaporator  divided  by  the  product  of  local  density  of  the 
working  fluid  vapor,  cross  sectional  area  of  the  vapor  space, 
and  the  local  heat  of  vaporization  of  the  working  fluid.  Due 
to  the  physical  boundary  at  the  end  of  the  pipe,  where  x  =  0, 
the  velocity  at  the  farthest  end  of  the  condenser  away  from  the 
evaporator  is  equal  to  zero.  The  boundary  conditions  for  the 
dependent  variable  of  velocity  are: 

when   X  =  0,   V  =  0.0   and 
when   X  =  L,   V  =  V 


e      p   A  h^ 
a   c  fg 
e      -^e 


Based  on  the  system  of  three  equations,  only  one  of  the  above 

boundary  conditions  is  required.   The  boundary  condition  at  the 

origin  is  used  for  the  iteration  scheme  and  V   is  used  to 
^  e 

normalize  velocity. 

The  temperature  and  pressure  at  the  farthest  end  of  the 

condenser  from  the  evaporator  are  known  from  laboratory  data. 

Assuming  thermodynamically  independent  substances,  the  total 

pressure  is  the  sum  of  partial  pressures  and  the  total  density 

is  the  sum  of  partial  densities.   The  saturation  pressure  of 

the  working  fluid  vapor,  P  ,  could  be  determined  from  tabulated 

a 

data  by  entering  the  saturation  pressure  table  with  condenser 
temperature.   An  example  of  this  table  can  be  found  in  Eubank 
[Ref.  8]  for  methanol  saturated  vapor.   The  noncondensible  gas 
pressure,  P,  ,  can  be  calculated  by  subtracting  the  working 
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fluid  vapor  pressure  from  the  total  pressure.   With  this  data, 
the  initial  mass  fractions  can  be  established  by  using  the 
following  procedure  from  Reynolds  and  Perkins  [Ref.  6]: 


X 


P 
a 


X 


total 


Pfa 


total 


A/i  Lbm  ^        .    , 

X  M    =   — r- = : — for  species   a 

'^a  a      Lbmole  mixture        ^ 


»/f  Lbm  ,-        .    , ,  , 

Xk^^k   =   TK^^i^  ^^..^,.v^   foi^  species  'b' 


'b  b      Lbmole  mixture 


CO 


X  w 
^a  a 


0      '^a  a   '^b  b 


03 


Xb"b 


''O      Xa"a  ■''<b"b 


At  the  end  of  the  condenser  closest  to  the  evaporator,  it  is 
assumed  that  the  concentration  of  noncondensible  gas  is  zero 
so  that  at  this  point  00   is  equal  to  one.   The  boundary  condi- 
tions for  the  dependent  variable  of  mass  fraction  for  the 
working  fluid  vapor  are: 


when  X  =  0 ,   w   =00   ,   and 

^0 


when  X  =  L,   oj   =1.0. 
a 
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C.   EQUATION  NONDIMENSIONALIZATION 

To  generalize  the  equations  and  facilitate  their  solution 
it  is  convenient  to  use  dimensionless  variables  defined  as 
follows . 

(a)  The  independent  variable,  condenser  length,  is  defined 
by : 

X 

"^   =   L 

where  x  is  the  distance  along  the  heat  pipe,  measured  from 
the  condenser  end  farthest  from  the  evaporator,  and  L  is  the 
total  condenser  length. 

(b)  A  dimensionless  temperature  can  be  defined  by: 

T  -T 


T  -T 
e    < 


where  T  is  the  local  temperature,  T  is  the  evaporator  tempera- 
ture, and  T  is  the  ambient  temperature.  In  this  analysis,  the 
first  and  second  derivative  of  temperature  with  respect  to  dis- 
tance are  required. 


dT   _   dT  de  dn   _   "^e   "^^de 


dx      d6  dn  dx        L    dn 


and 


d^T   ^   ^fdT-   ^   d_ /^e_^  d0_^drl 

^  2      dx^dx^      dn    L    dn  dx 
dx 


T  -T   -.2^ 
e   CO  d  6 


L    dn 
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(c)   The  second  dependent  variable,  working  fluid  mass 
fraction,  is  dimensionless  by  definition  and  its  numerical 
value  is  greater  than  or  equal  to  zero  and  less  than  or  equal 
to  one.   The  first  derivative  of  working  fluid  mass  fraction 
with  respect  to  distance  is: 


dco 

dco       -, 

1    dco 

a 

a   dn 

1         a 

dn 

dn    dx 

L    dn    * 

(d)   Velocity  of  the  working  fluid  vapor  is  made  dimension- 
less  by  normalizing  local  vapor  velocity  with  the  evaporator 
vapor  velocity,  V  .   The  nondimensional  velocity  and  the  first 
derivative  of  velocity  with  respect  to  distance  are  then 
defined  by: 


V 


V 


dV   ^   dV   dV*  di2   ^   ^  dV* 
dx      dV*  dn  dx       L   dn 


(e)   In  order  to  complete  nondimensionalizing  process,  the 
following  other  dimensionless  variables  are  used. 

(i)   The  working  fluid  vapor  density  is  normalized  with 
respect  to  the  evaporator  vapor  density. 


Sa*   =   Sa 


Sa 
e 
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(ii)   The  working  fluid  vapor  pressure  is  compared 
to  the  evaporator  vapor  pressure. 


T.  *      Fa 
Fa*   -   — 


Fa 
e 


(iii)   The  working  fluid  heat  of  vaporization  is  normalized 
with  respect  to  the  condenser  working  fluid  vapor  heat  of 
condensation . 


H 
H 


*   -   fg 


fg      H^ 
^       fg,c 


(iv)   The  working  fluid  vapor  viscosity  is  normalized 
with  respect  to  the  evaporator  vapor  viscosity. 


y*   =   -^ 
^e 


After  substitution  of  nondimensional  terms,  the  three 
equations  which  describe  operation  of  a  vertical  variable 
conductance  heat  pipe  become: 


(a)   -  ^(Sa^Jl^V^   d  D   d   Saj^  ^   ^   ^ 

dn   oj     dn      '  Lv  dn  co    d  n 
a  e     a     ' 


The  dimensionless  coefficient  from  Equation  (III. 8)  is 
D 


Nl  = 


LV  * 
e 
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(b)   ^  ^iJ-  +  ^-o (5-(Sa*(—  -1)  (e  + 


Sa^V  M,  V  a  e   °° 

e  e  be 


Lg  Sa^  _  d_  Sa*V*^>   ,  ,.*  d  ,Sa*V*. 

V    a  a  a 

e 


P   8  M   L 

— u*  V*  (III. 9) 

A   D  Sa   V   ^  vxxx.:?; 

c  V   e  e 


The  dimensionless  coefficients  from  Equation  (III. 9)  are 


Fa  R(T  -T  ) 

O  Q      OO 

N2   =   ^  ,     N3   =   ^.5 ,    N4   = 


2  '  2     '  T  -T  ' 

Sa  V  M,  V  e   ~ 

e  e  be 


J                                         P   8y   L 
N5   =   ^,   and   N6   = ^- 


,,2'  ^''^      ^'^  A  D  Sa  V   • 

V  c  V   e  e 

e 


-.2^  P  hL       H^     A  Sa  V  L 

(c)       ^— i  -  _s 0  _   fg/C  c   e  e   ^ 

^^    ,  2  K  A         K  A  (T   -T  )     fg 

dri         D  p   e    °° 


X  g-(      )   =   0  (III. 10) 

dn   co^ 

CI 


The  two  nondimensional  terms  of  interest  from  Equation  (III.  10 
are 


P   hL  H^     A   Sa   V   L 

,.T-7       s  J   ivTO       fq,c   c    e   e 

''^       =      -K^r-      ^^^   ^^   =    K  A  (T   -T  ) 

p  p   e    °° 


Numerical  evaluation  of  the  dimensionless  coefficients  is 
required  to  determine  the  relative  significance  of  each  term. 
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Two  experimental  runs  from  Batts  [Ref.  4]  are  used  to  evalu- 
ate the  dimensionless  coefficients.   A  low  powered  run 
(Q  =  22.9  watts)  with  methanol  as  the  working  fluid,  and  krypton 
as  the  noncondensible  gas,  is  chosen.   This  run  demonstrates 
an  example  of  the  heavier  molecular  weight  noncondensible  gas 
diffusing  through  a  large  portion  of  the  condenser.   A  second 
higher  powered  run  (Q  =  70  watts)  with  methanol  as  the  working 
fluid,  and  helium  as  the  noncondensible  gas,  is  chosen  because 
it  demonstrates  an  example  of  the  lighter  molecular  weight 
noncondensible  gas  occupying  a  small  portion  of  the  condenser. 
The  convection  heat  transfer  coefficient  for  the  heat  pipe  is 

Rti] 
set  equal  to  1.0  ^ ,  representing  free  convection.   A 

hr.ft.  °R 

R  t  u 
system  thermal  conductivity  of  approximately  4 . 56  ^ 

hr.ft.  °R 
results  from  a  parallel  circuit  combination  of  conductive 

thermal  resistances  for  the  stainless  steel  pipe  wall  and 

methanol  filled  stainless  steel  wire  mesh  wick.   Kane  [Ref.  5] 

calculated  representative  binary  mass  diffusion  coefficients 

2  2 

of  0.4  ft  /hr  for  krypton-methanol  and  1.9  ft  /hr  for  helium- 

methanol.   These  values  for  binary  mass  diffusion  are  also 

used  in  this  analysis.   The  experimental  values  that  are  used 

to  calculate  the  dimensionless  constants  are  listed  in  Table 

I.   The  values  of  the  dimensionless  constants  for  the  two 

experimental  runs  are  listed  in  Table  II. 

Based  on  the  tabulated  numerical  results.  Equation  (III. 9) 

is  modified  to: 
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TABLE  I 
NUMERICAL  DATA 


Constants  for  both  runs 


R  =  1545.33  ft  Lbf/Lbmole  °R 

2 
g  =  32.2  ft/sec 


Gas  constant 


Acceleration  of  gravity 


L  =  3.5  ft 


Length  of  condenser 


A   =  .0199  ft" 
c 


Cross  sectional  area  of 
vapor  space 


A   =  .001917  ft' 
P 


Axial  cross  sectional  area 
of  wick-wall  combination 


D   =  .15917  ft 

V 


Diameter  of  vapor  space 


P   =  .524  ft 
s 


Periphery  of  outside  of 
pipe 


P   =  .500  ft 

V 


Periphery  of  vapor  space 


h  =  1.0  Btu/hr  ft  °R 


K  =  4.5585  Btu/hr  f t °R 


Convection  heat  transfer 
coefficient 

Thermal  conductivity 


Low  powered  krypton/methanol 


Q  =  22.9  watts  =  78.1348  Btu/hr 


T   =  592°R 
e 


T   =  543°R 
c 


T   =  536°R 


D  =  0.4  ft  /hr 


M,  =  8  3.8  Lbm/Lbmole 
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TABLE  I  (CONTINUED) 


V   =  .041339  ft/sec 
e  ' 

Fa   =  10.24  Lbf/in^ 


H^     =  500.21  Btu/Lbm 
fg,c 

Sa   =  0.054871  Lbm/ft^ 

y   =  224.502  x lo"^  Lbf  sec/ft^ 
e 


High  powered  helium/methanol 

Q  =  70  watts  =  238.84  Btu/hr 


T   =  640°R 
e 


T   =  540°R 
c 


T   =  536°R 


D  =  1.9  ft^/hr 


M,  =  4.003  Lbm/Lbmole 


V   =  .050485  ft/sec 
e 

Fa   =  28.52  Lbf/in^ 
e 


H^     =  501.11  Btu/Lbm 
fg,c 

Sa   =  0.14446  Lbm/ft^ 

y   =  243.993  x 10~^  Lbf  sec/ft' 
e 
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TABLE  II 
DIMENSIONLESS  CONSTANTS 

Low  Powered  Krypton/Methanol        High  Powered  Helium/Methanol 


Nl   =  j^  =    7.679  X  10  ^ 

1j  V 


Nl  =  j^   =  2.9869  X  10  -^ 


N2   = 


Fa  f-  Fa 
%-  =      506.355x10       N2   =   —   ® 


Sa      V 
e      e 


Sa   V^ 
e    e 


359.166  X 10 


R(T    -T    )                                            ^                           R(T    -T    )  ^ 

N3      =      ^    ^°°         =      19.4581x10  N3      =      ^-^ =       507.224x10^ 

be  be 


N4       = 


T    -T 
e      ^ 


=      9.571 


N4      = 


T    -T 

e      -^ 


=       5.154 


N5      =      ^     =      65948.4 
V 


N5      = 


-      '^      =      44218.0 


V 


N6       = 


P      8  y       L 
V        e 

A     D    Sa      V 
c     V       e     e 


=      14.086  N6 


P      8  VI      L 
V        e 

A     D    Sa      V 
c     V       e     e 


4.76143 


P     hL 

N7      =      ^— =      734.553 

K    A 


P      hL 
N7      =      -^ =      734.553 


H^  A      Sa      V      L 

MR      -         fg/^      c         e      e 
^"      ~  K   A    (T    -T    ) 

p      e      °° 


N8       = 


H^  A      Sa      V      L 

fg,c      c         e      e 

K   A    (T    -T    ) 
p      e      °° 


=      581.362 


=      1008.30 


Fa  ,^  *    R(T  -T  )  ,        ,            T 

e  dFa*  .     e   <»'  .d  ,_,  ^  ,  1   i  \  /a  ,    °° 

-— 2  -d^  -^    ,,   2   t^^SaM—  -1)(9  4-^r^^ 
Sa  V           M,  V      '       a          e   = 


-)] 


e  e 


be 


Lg  Sa' 


V 


LO 


(III. 11) 


The  viscous  term  and  the  two  inertia  terms  are  neglected  be- 

1  3        fi 

cause  their  order  of  magnitude  is  10   compared  to  10   and  10 

for  the  remaining  terms.   The  nondimensional  equations  (III. 8), 

(III. 10),  and  (III. 11)  are  used  to  develop  the  analytical  model 

of  the  vertically  oriented, variable  conductance  heat  pipe. 

D.   FORMULATION  OF  THE  NUMERICAL  MODEL 

Equation  (III. 8)  is  integrated  and  is  solved  for  the  first 
derivative  of  the  working  fluid  mass  fraction.   Integrating 
both  sides  of  Equation  (III. 8)  gives: 


Sa*V*  _^  g^^^^  _  ^-L  Sa* 


dto 


0) 


CO 


dn 


(III. 12) 


The  lower  endpoint  is  used  (n  =  0)  to  evaluate  the  constant 
of  integration.   When  n  =  0,  V*  at  the  end  of  the  pipe  equals 
zero.   The  constant  of  integration  then  becomes: 


C   =   -Nl 


'0 


doj 


Sa   dn 


(III. 13) 


By  replacing  the  value  for  the  constant  of  integration  and 

solving  for  dw  /dn ,  the  final  version  of  Equation  (III. 12)  is 
a 
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dco 


-V* 


dn  Nl  Nl 


V*a)  CO         p^      dw 

a  a         0  a 


Sa*    Sa        dn 


(III. 14) 


The  derivative  in  Equation  (III. 11)  is  expanded  to  get  the 
following . 


|^(SaMji-l)(9  +N4)) 

3. 


^ (—  -1)  (e  +N4)  +  Sa*  ^  -AiQ   +N4) 

dn   CO  2   dn 


CO 


+  sa*(^-l)  ^ 

CO       dn 


(III. 15) 


The  chain  rule  is  applied  to 


dSa* 

dn  • 


dSa*   ^   dSa*  dSa  dT  d9 
dn       dSa   dT  dG  dn 


T  -T 
,  e   °°>dSa  de- 

^  Sa   ^  dT  dn 
e 


(III. 16) 


dFa- 


Applying  the  same  procedure  to  — -= ,  the  result  is: 


dFa*   ^   /^ej^>  dPa  d9 

dn       ^  Fa   ^   dT  dn 
e 


(III. 17) 


d  fl 
By  combining  terms  and  solving  for  -r—,    the  final  form  for 

equation  (III. 11)  is: 
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(a)   The  coefficient  of  -r—   is  set  eaual  to  G  (w  ,6). 

dn  a'  ' 


N2 (T  -T  )   .^     N3 (T  -T  )   ,^    , 
G   =   ^ ^  +  ^rf ^(—  -1)  (6  +N4) 


Fa       dT 
e 


Sa       dT  'co 
e  a 


+  N3(Sa*)  {-^  -1)  . 


OJ 


(III. 18) 


d  ft 
(b)   Solving  for  3—,  the  result  is 

^     dn 


de   _   rN5Sa*  _^  N3Sa*  ^'^a  ._  _^^^..  ,  .^  .    .. 

3 —   =   I +  T —  —5 — (6  +N4)J/G(uj  ,0) 

dn         w         2    dn  a ' 

a      w 


(III. 19) 


Equation  (III. 10)  is  integrated  and  the  result  is: 


|i-N7   /   edn  -  N8  H*  §i^  =  c  . 

dn      f.^  fg   CO 


(III. 20) 


Using  the  boundary  condition  of  n  =  0  where  V*  =  0.0,  the 
constant  of  integration  is: 


C   = 


d^ 

dn 


(III. 21) 


The  final  version  of  Equation  (III. 10)  which  is  solved  for  V* 
is  : 


V*   =   (^  w  -N7(/edn)a3  -03   ^^ 


duu 
n  a 


a    a  dn 


0 


)/(N8H*  Sa*)  .  (III. 22) 
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The  final  system  of  three  equations  with  nondimensional 
variables  is: 


do)  dco 

(a)       Equation    (III. 14),    -^(v*,oj    ,-^ 

dn  a     dn 


) 


dw 


T7+         V*co  0)         p^      dco 

-^^*  a    ,       a      ^0  a 


■V^ 


dn 


Nl  Nl  S^*    Sa        dn 

Q  e 


(b)       Equation    (III. 19),    J^(co    ,6,-5-^); 

dn       a  dn 


de_ 
dn 


rNSSa*  ^N3Sa*    ^^a  ,.  ,^,.>  T   ,^  ,         „, 

+ T ^-r   9+N4)    /G    CO    ,9) 

CO  2  dn  a 

a  CO 


where 


N2 (T    -T    )      ,^  N3(T    -T    )     ,^         , 

n  /         a\  e      °°      dFa    ^  ^    e      °°      dSa  ,1       1  a  /a  ,>t/i  ^ 

^(^a'Q^   =   Fi ^T  -^  Si -df(—  -.^^^Q-'N^) 


+  N3(Sa*)  i-!^  -1) 


CO 


(c)   Equation  (III. 22),  V*  (co  ,  f  Odn,^,^ 
^  a  ^^      'dn  dn 


V*   =   (^00  -N7(fedn)co   -co   ^^ 


dn   a 


a    a  dn 


0 


)/(N8H|gSa*)  . 


The  boundary  conditions  for  this  set  of  equations  are  the 
following : 
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(a)   For  nondimensional  temperature 


n   =   0.0 


T  -T 
c    <^ 

T   -T 
e    <: 


n   =   1.0 


=   1.0 


de_ 
dn 


0 


=   a 


The  constant  a  is  not  known  but,  based  on  laboratory  data,  can 
be  estimated  and  through  an  iteration  process  the  correct  value 
can  be  obtained.   A  proposed  iteration  scheme  is  presented  in 
Chapter  IV. 

(b)   For  working  fluid  vapor  mass  fraction: 


n   =   0.0 


00    =   to 
a      a, 


n  =  1.0 


dco 


dn 


=   3 


0 


OJ 


=   1.0 


The  constant  3  is  not  known,  but  can  also  be  estimated  for  an 
initial  value  and,  through  an  iteration  scheme,  a  correct  value 
can  be  obtained. 

(c)   For  working  fluid  vapor  dimensionless  velocity: 


n   =   0.0 


V 


*   — 


0.0 


43 


E.   AUXILIARY  EQUATIONS 

A  means  to  derive  saturated  working  fluid  heat  of  vapori- 
zation, density,  and  pressure  is  required  before  a  numerical 
solution  is  attempted.   Using  methanol  as  the  working  fluid  and 
data  from  Eubank  [Ref.  8]  for  saturation  conditions,  the 
following  functions  are  derived. 

1 .  Polynomial  for  the  Working  Fluid  Heat  of  Vaporization 
Tabulated  values  of  the  heat  of  vaporization  for  methanol 

are  plotted  on  Figure  3.1  with  a  solid  line.   A  Fortran  program, 
polyfit,  is  used  to  obtain  the  coefficients  of  the  following 
polynomial : 

H^  (T)   =   aT^  +  bT  +  c  (III. 23) 

fg 

where  a  =  -9.96  x lo"^  Btu/Lbm ( °R) ^ ,  b  =  .730  Btu/Lbm(°R),  and 
c  =  397.99  Btu/Lbm.   The  values  of  H^  (T)  using  this  polynomial 
are  plotted  on  Figure  3.1  with  a  dotted  line. 

2 .  Polynomial  for  the  Working  Fluid  Vapor  Saturation 
Density 

Tabulated  values  for  methanol  saturated  density  are 

plotted  on  Figure  3.2  with  a  solid  line.   Again  the  Fortran 

program,  polyfit,  is  used  with  the  following  modification: 


'*   =  7^r^     where   T   ^   =   SOO'^R  and   C   =   e 
T   ^  ref 

ref 


A  polynomial  in  the  variable  C  is  used  to  describe  the 
saturation  density. 
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The  polynomial,  Sa(5),  is  then: 


Sa(C)   =   aC"^  +  b^^  +  cC  +  d  (III. 24) 


where  a  =  .1293  Lbm/ft^,  b  =  -1.0475  Lbm/ft^,  and 

c  =  2.8868  Lbm/ft^  and  d  =  -2.7037  Lbm/ft^.   The  value  of 

dSa/dT  is  also  required  for  the  numerical  solution.   This  is: 


^  =  (3a^^  +  2bC^  +  cU/T^gf  (III. 25) 


3 .   Polynomial  for  the  Working  Fluid  Vapor  Saturation 
Pressure 

Tabulated  values  for  methanol  saturation  pressure  are 
plotted  on  Figure  3.3  with  a  solid  line.   Polyfit  is  used  with 
the  same  modification. 

The  polynomial  Fa(^)  is: 


Fa(C)   =   aC^  +  b?^  +  c?  +  d 


where  a  =  24.5994  Lbf /in^ ,  b  =  -191.2472  Lbf/in^, 
c=  498.4940  Lbf/in^,  and  d  =  -439.9883  Lbf/in^.   The  first 
derivative  of  Fa  (5)  with  respect  to  temperature  is  equal 
to : 


^  =  (3a?^  +  2bC^  +  ^^^/^ref  '  (III. 26) 
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IV.   DISCUSSION  AND  RECOMMENDATIONS 

A.   DISCUSSION 

The  final  version  of  the  model  for  a  variable  conductance 
heat  pipe  with  a  vertical  orientation  is  a  system  of  three 
ordinary,  non-linear,  coupled,  differential  equations.   These 
three  equations  are  (III. 14),  (III. 19),  and  (III. 22). 

The  remaining  challenge  for  the  numerical  model  is  to  pro- 
pose an  iteration  scheme  that  can  solve  this  system  of  three  equa- 
tions.  A  special  thank  you  goes  to  the  collaboration  of  Dr.  M.D. 
Kelleher  and  Dr.  K.T.  Yang  for  the  following  iteration  technique. 

The  steps  of  the  iteration  procedure  are  as  follows: 

(a)   To  start  the  iteration  process,  assume  a  form  for  the 

relationship  between  6  and  co   such  that  co   =0.6+0^.   This 
^  a  a     1      2 

linear  relationship  is  reasonable  because: 

(i)     When  n  =  0 ,  both  co   and  9  are  values  close  to 

a 

zero . 

(ii)    When  n  =  L  both  co   and  9  are  equal  to  1 . 0 . 

3. 

(iii)   This  step  is  required  for  the  first  iteration 
only. 
The  equations  for  this  step  are  then: 


00    =   C-^e  +  C2  (IV. 1) 


dco  -,„ 

-^   =   C,  ^  .  (IV. 2) 

dn      1  dn 
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d  ft 
Then,  using  Equation  (III. 19),  solve  for  ^.      This  establishes 

a  first  degree  equation,  g— ,  in  terms  of  6  only  which  can  then 

be  integrated  numerically  to  establish  a  first  iteration 

relationship  for  G(ri). 

dfl 

(b)  With  this  relationship  for  9  and  ^—   defined  from  step 

^  dn  ^ 

(a)  above.  Equation  (III. 22)  can  be  solved  for  V*. 

(c)  Once  step  (b)  is  completed,  Equation  (III. 14)  can  be 

solved  for  dco  /dn .   The  integral  of  dco  /dn  is  then  used  to 
a  a 

define  w  .   The  process  can  then  be  repeated. 

In  summary,  the  first  iteration  goes  through  steps  a-c 

and  the  remaining  iterations  go  through  steps  b-c  until  there 

are  no  further  changes  in  the  form  of  9  ,  oj  ,  and  V*.   In 

a 

theory,  the  outcome  of  this  numerical  iteration  is  a  table 

which  gives  dimensionless  temperature,  9;  velocity,  V*;  and 

mass  fraction  of  the  working  fluid  vapor,  co  ,  versus  non- 

a 

dimensional  condenser  length,  r\ . 

B.   RECOMMENDATIONS 

1.  This  numerical  analysis  was  written  to  yield  results  via 
a  computer  program.   The  Continuous  System  Modeling 
Program  III  [Ref.  9]  is  a  computer  package  that  could 

be  used  to  solve  the  system  of  equations.   The  CSMP 
[Ref.  9]  includes  a  powerful  integral  statement  which 
carries  out  a  Fourth  Order  Runge-Kutta  Method. 

2.  This  model  assumes  constant  values  for  the  binary  mass 
diffusion  coefficient  and  the  convection  heat  transfer 
coefficient.   Once  numerical  values  are  obtained  from 
the  iteration  scheme,  the  accuracy  can  be  improved  by 
using  analytical  expressions  instead  of  constants  for 
the  binary  mass  diffusion  coefficient  [Ref.  10]  and  the 
convection  heat  transfer  coefficient  [Ref.  11]. 

3.  An  attempt  could  be  made  to  model  the  variable  conduc- 
tance heat  pipe  that  is  not  in  a  vertical  orientation. 
Kelleher  [Ref.  3]  shows  that  a  stratification  of  the 
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working  fluid  and  noncondensible  gas  occurs  when  the 
heat  pipe  is  in  the  horizontal  position.   This  position 
makes  the  transport  processes  inside  the  pipe  three 
dimensional,  which  would  require  a  finite  difference 
computer  model  for  solution. 
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